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We study theoretically the current-induced magnetic domain wall motion in a metallic nanowire 
with perpendicular magnetic anisotropy. The anisotropy can reduce the critical current density of 



*^J . the domain wall motion. We explain the reduction mechanism and identify the maximal reduction 

conditions. This result facilitates both fundamental studies and device applications of the current- 



induced domain wall motion. 
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Spin-polarized electrical currents in ferromagnets can transfer their spin angular momen- 
tum to local magnetizations via the s-d exchange interaction and generate torques [l|, y] 



) received considerable attention in 
7| and applications [a |9|] . 



on local magnetizations. This spin transfer torque (S r . 
view of both fundamental physics research 3|, U, |5|, |g, 

In a ferromagnetic nanowire, the STT can generate motion of magnetic domain walls 
(DWs). For conventional metallic ferromagnetic nanowires, which have the in-plane mag- 



netic anisotropy (IMA), experiments 10 



11 



121 ] found such current-induced DW motion 



when the current density J in the nanowire is larger than a certain threshold value J c of 
the order of 10 8 A/cm 2 . This value is too hi gh; At such current densities, the Joule heating 



13 



14 



15| , making fundamental studies of the 



generates considerable thermal fluctuations _, ., ...... 

STT difficult. Furthermore device applications [16| require J c < 10 7 A/cm 2 at room tem- 
perature. Thus both for fundamental studies and device applications, it is crucial to reduce 



Recently there are experimental 13] and theoretical 18J indications that J c may be con- 
siderably lower in a metallic nanowire with the perpendicular magnetic anisotropy (PMA). 
However it remains unclear how the PMA can lower J c . We aim to answer this question in 
this Letter. 

We consider a nanowire with the wire width w along the y-axis and thickness t along the 
z-axis (w > t). We use the Landau-Lifshitz-Gilbert (LLG) equation with the STT term, 



rh = —7m x H e flf + am x m + bj(J ■ V)m 
— cjm x (J • V)m, 



(1) 



where m is the unit vector of the local magnetization, 7 is the gyromagnetic ratio, a is the 
Gilbert damping parameter, J is the unit vector of the local current density, and H e g is 
the effective magnetic field, bj = P^bJ '/eM s is the magnitude of the adiabatic STT (5(, 
where e is the electron charge, P is the spin-polarization of the ferromagnet, \xb is the Bohr 
magneton, and M s is the saturation magnetization, cj is the magnitude of the nonadiabatic 
STT |6|, lZI] with the non-adiabaticity represented by the dimensionless parameter (3 = cj/bj. 
(3 is independent of J and estimated [6| to be of the order of 1CT 2 . 

To get an insight into the main physics of the PMA, we first develop an analytical 
model based on a one-dimensional (ID) approximation. Its results will be later verified by 
performing the micromagnetic simulations of the LLG equation [Eq. (CQ)], which are known 



to provide a reliable description of nanoscale magnetization dynamics 19|, [20|. 
For a ferromagnet with the PMA we have 

H - m^ + ^r z + dipole ' (2) 

where A is the exchange stiffness constant and Ky is the PMA constant that allows the easy 
axis (along the z-axis) to be perpendicular to the wire-plane (x-y plane). To describe the 
demagnetization effects, we consider the magnetostatic dipole-dipole interaction field given 
by Hdipoie(r) = M g J d 3 r'N(r — r')m(r'), where the components of the matrix iV are given 
by N xx (r) = — [1 — 3x 2 /|r| 2 ]/|r| 3 , N xy (r) = 3xy/\r\ 5 . Other components of N are defined in 
a similar way. 

We also assume that the DW maintains the following shape during the DW motion; 
m z (r, t) = tanh[(s — q)/X], m x (r,t) = cos^sech[(x — q)/X], m y (r,t) = sin^sech[(x — q)/X], 
where A is the equilibrium DW width obtained from ID micromagnetic simulations. In 
this rigid DW motion approximation |2l| the DW dynamics is described by two dynamical 
variables, the DW position q(t) and the DW tilting angle tj){t). 

By using the procedure developed by Thiele 22|, one can then derive, from the LLG 
equation, the equations of motion for the two collective coordinates q and ip, 

X^-aq = Cj- ( 7 A/2M s )/ pin , (3) 

q + aX^J = -bj - ( 7 A/M s )K d sin 2^, (4) 

where the pinning force / pin is related to the DW energy per unit cross-sectional area u tot 
(/pm = —du tot /dq) representing pinning potential in the presence of extrinsic defects in a 
nanowire. Here K& is the effective wall anisotropy given by 

K d = K y - K x , (5) 

where K-i = —^^ JJ d 3 rd 3 r'Na(r — r')sech|sechy (i = x,y,z) and S is the cross-sectional 
area. K& represents the magnetostatic energy difference between two types of transverse 
DWs, the Bloch DW (m || e y at the DW center) and the Neel DW (m || e x at the DW 
center) . 

Before we demonstrate its implications for a general case, we first consider a defect-free 
nanowire (/ pin = 0) with cj = 0. J c in this case is given [5|] by J" 



Till 

' (: 



J? , $±\U. (6) 
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Figure Q^a) shows w and t dependence of J™. Note that J c falls below 10 7 A/cm 2 in a 
wide range of w and t. Since all material parameters used in Fig. [1] are similar to those for 
permalloy except for the PMA constant K V) this reduction in J™ should be attributed to the 
PMA. To check the validity of this prediction, we also perform micromagnetic simulations 
of the LLG equation and excellent agreement is found [Fig. [TJb) upper panel] . 

This reduction in J™ becomes especially effective when w is tuned to a ^-dependent 
special value w*(t), at which K d reverses its sign [Fig. QJb) lower panel] and near w*(t), J™ 
(oc \K d \) is strongly suppressed. The sign reversal of K d implies that w* is the equilibrium 
phase boundary between the Bloch DW and Neel DW. For transverse DWs in an IMA 
nanowire, on the other hand, K d is given by K z — K y and since K z is always larger than K y , 
K d in the IMA case is always positive in a conventional nanowire geometry with w > t. This 
difference between a PMA nanowire and an IMA nanowire illustrates a crucial role played 
by the PMA. 

Next we consider a general case with / P j n 7^ and cj 7^ 0. After some calculation, one 



can obtain an upper bound J" p of J c [231 ]. 



J c up ee max [min(4 n , J c ex ), a/3J c ex ] , (7) 

where J^ x = ( r y\e/2Pfi B )(f™*/j3) and f^ x represents the maximum value of / pin . The 
dashed line in Fig. [T](c) shows J" p as a function of f^ x for a PMA nanowire. For the 
case J™ = 1.6 x 10 6 A/cm 2 , Fig. [T]^c) also shows J c determined from numerical simulations 
of Eqs. ([3]) and (jlj) with the pinning potential energy w tot modelled by a finite ranged 
harmonic potential [Fig. [T]^c) inset]. A few remarks are in order. Firstly, both J c and J" p 
exhibit plateaus near J™ in a wide range of f^ x . Secondly, J c depends on f3 only in the 
weak pinning regime (/™^/2M s < 1 Oe) and the (3 dependence essentially disappears in 
the intermediate (plateau) and strong (above plateau) pinning regimes. This behavior is 



consistent with the prediction of Eq. ([7]). Thirdly, a recent experiment [24| with a PMA 
nanowire finds the depinning magnetic field of about 500 Oe for a field-driven DW motion. 
When this value is used as an estimation of /"j^ x /2M s , one finds J" p ~ J™ ~ 10 6 A/cm 2 . 
Thus Fig. [T^c) demonstrates that the reduction of J™ via the PMA indeed leads to the 
reduction of J c . As a comparison, results for an IMA nanowire are also given in Fig. QJc). 
Differences from the PMA case are evident. 

Next we present micromagnetic simulation results of the LLG equation. Various sources 



of /pin are examined. Figure^a) shows J c obtained from the ID LLG equation for a situation 
where the magnitude of the PMA constant K\j fluctuates from its bulk average value Ku,o 
with the maximum deviation given by Vo- Note that the result is remarkably similar to that 
in Fig. QJc). In good quality PMA samples, Vq/Kjjq is reported [25] to be less than 0.1, for 
which we obtain J c rj 10 6 A/cm 2 . Figure MJo) shows the effect of a notch investigated with 
the two-dimensional (2D) LLG equation. In a wide range of w, J c falls below 10 7 A/cm 2 
despite the notch formation. Note that for w > 80 nm, J c decreases as the notch depth 
Sw increases. This strange behavior is not due to the locally enhanced current density near 
the notch, since this effect should be stronger for w < 80 nm. Instead it is due to the fact 
that J™ is determined by an effective wire width that a DW senses. When J c is plotted as a 
function of w — 5w/2, an estimation of the effective width, this strange behavior disappears 
and J c is now almost independent of 8w, in agreement with the prediction J" p = J™ in 
the plateau range in Fig. QJc). Figure [2]^c) shows J c for a PMA nanowire with the edge 
roughness and with the PMA fluctuations. Although values of J c are somewhat scattered 
with the realizations of the randomness, J c still remains below 10 7 A/cm 2 in a wide range 



of the average width w ave [26[] . 



Here we remark that all demonstrations for the reduction of J c assume the proper tuning 



of w and t to achieve the reduced J™. A recent experiment on the PMA nanowire [24| found 
J c — 1.0 x 10 8 A/cm 2 without such tuning. We suggest that the tuning of w and t can reduce 
J c . Another experiment |l7|] found indications of the enhanced STT efficiency in a PMA spin 
valve. However the measurement was still restricted to the thermally assisted creep regime 
with extremely low DW velocity (average t>d„ < 10 -8 m/s). According to our calculation 
(not shown), much higher velocity (t>d w ~ 10 m/s) can be achieved at J ~ 10 7 A/cm 2 



if w and t are properly tuned. Finally the report 27| of the reduced J c in ferromagnetic 
semiconductors is yet limited to low temperatures (~ 100 K) while the reduction scheme 
presented in this Letter does not require low temperature operation. 

In summary, we have clarified the mechanism by which the PMA can drastically reduce 
J c . The geometrical tuning is important to maximize the reduction by the PMA. When 
properly tuned, the dependence of J c on (3 and the DW pinning force f pin is very weak. 
This result solves the large thermal fluctuation problem and also makes feasible nanoscale 
magnetoelectronic devices J8j, |9( based on the current-induced DW motion. 
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FIG. 1: (Color online) (a) J™ from Eq. Q as a function of w and t. (b) Upper panel: J™ defined 
in Eq. Q vs J c obtained from micromagnetic simulations of the ID and 2D LLG equation. Lower 
panel: K&/M^ as a function of w. t = 10 nm in both panels, (c) The dashed (dash-dotted) 
line shows J c up [Eq. ©] for a PMA (IMA) nanowire with (5 = a and J™ = 1.6 x 10 6 (3.15 x 10 9 ) 
A/cm 2 . Symbols shows J c obtained from numerical simulations of Eqs. ([3]) and (j3|). The result is in 
reasonable agreement with J" p for (5 = a. Inset: Spatial profile of Utot with qo = 3 A. The following 
parameters are used: a = 0.02, and P = 0.7, A = 1.3 x 10~ 6 erg/cm, #"[/ = 1.5 x 10 6 (0 x 10 6 ) 
erg/cm 3 , and M s = 400(800) emu/cm 3 for PMA (IMA) nanowire. 
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FIG. 2: (Color online) Micromagnetic simulation results of the LLG equation, (a) Effects of the 
magnitude fluctuation of the PMA constant K\j for w = 77 nm and t = 10 nm. Upper inset: 
Spatial profile of Kjj with qo = 37.5 nm. (b) Effects of a notch. Inset: Schematic of a notch. 
(c) Combined effects of edge roughness (2.5 nm for each edge) and PMA fluctuations (Gaussian 
magnitude fluctuations of 5% and direction fluctuations of 5° for each cell of size 2.5 nm x 2.5 nm 
x t nm). For each wj ave , three realizations of the randomness are considered (The dotted line is a 
guideline). (3 = a and t = 10 nm. All other parameters are the same as those in Fig. [TJ 



10 



